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00 . Abstract 

o 

We study corrections suppressed by one power of the soft gluon energy to the 
resummation of threshold loErarithms for the Dreh-Yan cross section and for 
>< : Deep Inelastic structure functions. While no general factorization theorem 

j^ ■ is known for these next-to-eikonal (NE) corrections, it is conjectured that at 

least a subset will exponentiate, along with the logarithms arising at leading 
power. Here we develop some general tools to study NE logarithms, and we 
construct an ansatz for threshold resummation that includes various sources of 
NE corrections, implementing in this context the improved collinear evolution 
recently proposed by Dokshitzer, Marchesini and Salam (DMS). We compare 
our ansatz to existing exact results at two and three loops, finding evidence for 
the exponentiation of leading NE logarithms and confirming the predictivity 
of DMS evolution. 
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1 Introduction 

Sudakov resuminations are established in perturbative QCD for all logarithmic con- 
tributions, to leading power in the total momentum fraction carried by soft gluons. 
To illustrate this fact, consider as an example threshold resummation for the Drell- 
Yan process, or for a similar electroweak annihilation cross section at the hard scale 
Q. In this case, large logarithms arise in the hard partonic cross section when the 
total available center-of-mass energy, s, is only slightly larger than the mass Q^ 
of the selected electroweak final state. Gluon radiation into the final state is then 
forced to be soft, as gluons carry (at most) a total energy (1 — z)s, with z = Q"^ /s. 
As a consequence, perturbative contributions at order a" are enhanced by large 
logarithms in the form of 'plus' distributions, up to [ln^"~^(l — z)/{l — z)\ . Upon 
taking a Mellin transform, these distributions turn into powers of logarithms of 
the Mellin variable A^, conjugate to 2, up to In " A^. All these contributions can 
be resummed [H E], and they display a nontrivial pattern of exponentiation: the 
logarithm of the cross section in Mellin space, in fact, is enhanced only by single 
logarithms, up to In""*"^ A^ at order a". 

It has been understood since the early days of QCD [3j| that at least some non- 
logarithmic contributions (terms independent of A^, which are Mellin conjugate to 
virtual corrections proportional to 5(1 — z)) also exponentiate. In fact, Ref. [1] later 
proved, at least for electroweak annihilation and DIS, that all such contributions 
can be organized in exponential form. One may naturally wonder to what extent 
this pattern of exponentiation can be extended beyond leading power in the Mellin 
variable A^, or in the soft gluon energy fraction 1 — z. 

There are several problems in attempting to extend the resummation formalism 
beyond leading power in A^, or 1 — z. Indeed, resummation can be understood to 
be a consequence of Sudakov factorization, as discussed in [5]. To leading power 
in A^, it can be shown that the Mellin moments of the cross section factorize into 
distinct functions responsible for infrared and coUinear enhancements, times a hard 
remainder which is free of logarithms. Exponentiation follows from evolution equa- 
tions that are dictated by this factorization. To date, no proof of such a Sudakov 
factorization is available beyond leading power in A^. Part of this issue is the fact 
that, in order the achieve exponentiation, the phase space specific to the observable 
at hand, in the threshold limit, must itself factorize; this is achieved at leading power 
by taking the Mellin transform, thanks to the fact that the observable (essentially 
1 — z for the inclusive Drell-Yan cross section) is linear in soft gluon energies to 
leading power in 1 — z. Again, this simple property is lost beyond leading power. 

Not withstanding these obstacles, there is intriguing, if scattered, evidence that 
some of the mechanisms that lead to the resummation of leading power logarithms 
are still operating at next-to-leading power. Theoretically, evidence in this direction 
is provided by the Low-Burnett-Kroll theorem [01 [7|, which states that (in QED) 
cross sections involving soft photons can be expressed in terms of radiation-less 
amplitudes not only at leading power in soft gluon energies (which corresponds to 



the breinsstrahlung spectrum and to the eikonal approximation), but also at next- 
to-leading power. For such cross sections radiation is simply related to classical 
fields, and one expects some form of soft photon exponentiation to hold. In QCD, 
direct application of Low's theorem is complicated by the presence of coUinear di- 
vergences |8j|, but one may still expect it to be relevant for soft emissions. 

At a more practical level, one may observe that resummed cross sections are 
expressed in terms of integrals of certain anomalous dimensions, with integration 
limits dictated by the phase space available for soft radiation, and with the running 
coupling evaluated at the typical transverse momentum of the first gluon emission. 
These kinematical quantities are evaluated in the threshold limit, and one may 
expect that correcting their values in order to make them accurate at next-to-leading 
power in the soft momentum should lead to a physically meaningful improvement 
of the resummation. 

This kind of reasoning has led to attempts to include certain sub-eikonal effects 
in practical implementations of Sudakov resummations, mostly in view of gauging 
the theoretical uncertainty of the resummation [9]. Typically, this involves including 
subleading terms in the coUinear evolution kernel into the resummation, which is 
particularly appealing for Drell-Yan and related cross sections, where the entire 
singularity structure is determined by initial state soft and coUinear radiation. This 
was applied in the case of Higgs production in Refs. P, [HI [TOl [I2], and for prompt 
photon production in Ref. [13]. 

More recently, following the evaluation of coUinear evolution kernels at three 
loops [H], a bold suggestion has been put forward by Dokshitzer, Marchesini and 
Salam (DMS) pLSj, who proposed a modified evolution equation for parton distri- 
butions, based on the idea that the proper ordering variable in the coUinear shower 
should be the lifetime of parton fiuctuations rather than the gluon transverse mo- 
mentum. This modified evolution has remarkable consequences: it explains a previ- 
ously mysterious numerical coincidence observed by [H] , and it connects eikonal and 
sub-eikonal terms in the splitting function in a nontrivial way, consistent with the 
idea that all evolution effects which are non- vanishing as 2; — > 1 should be determined 
at one loop, with an appropriate definiton of the coupling. The DMS proposal has 
later been refined by Basso and Korchemsky [16], who traced the recursive relation 
which determines the coUinear anomalous dimension to the conformal invariance of 
the classical theory, and its breaking by the (3 function. The relations connecting 
eikonal and next-to-eikonal terms for parton evolution are then generalized to higher 
twist operators as well. 

In this note, we begin to develop a systematic approach for the inclusion of next- 
to-eikonal terms in the resummation, inspired by the results of [15] and by the earlier 
work of [9]. We begin, in Section [2], by briefiy reviewing the DMS approach, and 
describing how we intend to implement it in the context of Sudakov resummation. 
There, we also introduce some simple tools and definitions to evaluate the integrals 
that appear in resummed exponents to the desired accuracy. Then, in Section [31 we 
propose an ansatz to include in the resummation all next-to-eikonal effects that can 



be argued to be under theoretical control. We do this for the Drell-Yan cross section 
and for the Deep Inelastic structure function F2. It is clear from the outset that our 
ansatz controls only a subset of all next-to-eikonal terms in the cross section: indeed, 
it may well be that not all such terms can be organized in exponential form. We be- 
lieve however that the terms we include are physically well motivated, so we expect 
our ansatz to reproduce with reasonable accuracy higher order perturbative results, 
based on the evaluation of the exponent at lower orders. We proceed to test this 
expectation by comparing the results of expanding our proposed resummed expres- 
sions with the known exact results at two loops for the Drell-Yan cross section [17] , 
and at two and three loops for DIS [HI [19] . 

As we will outline in our discussion, in Section HJ the results of this comparison 
are consistent with the assumption that at least leading next-to-eikonal logarithms 
do exponentiate, for all color structures. Furthermore, the implementation of the 
DMS approach reproduces with considerable accuracy (though not exactly) certain 
classes of subleading next-to-eikonal logarithms which could not have been gener- 
ated by the standard resummation. We believe that these results are encouraging 
regarding the possibility that next-to-eikonal logarithms could be understood and 
organized to all orders, an effort which will ultimately require a full analysis of soft 
gluon effects beyond the eikonal approximation. 

2 Tools for next-to-eikonal resummation 

The task of probing the extension of the resummation formalism beyond the eikonal 
approximation requires both conceptual and practical tools. In this section we 
describe briefly the main conceptual progress that we are going to employ, which is 
the idea, put forward by DMS, that all NE terms in collinear evolution trace their 
origin to one loop effects, phase space, and the choice of an appropriate, physically 
motivated coupling. We present the DMS equation, and we show how it can be 
solved in exponential form, just like ordinary collinear evolution, to NE accuracy. 
Next, making use of a technique developed in [2U], which we generalize to NE level, 
we present some simple results for the generic integrals that may appear in NE 
resummed cross section to any perturbative order. 

2.1 The DMS evolution equation and its solution 

Consider first the familiar collinear evolution equation for the non-singlet quark 
density 

d f dz / X \ 

/^^^^(^'/^^) = / y^i^'^y Pqq{z,as{^i'^)) ■ (1) 



As is well known, this simple convolution can be turned into a product by taking a 
Mellin transform, 

/x^ A q^N, /i^) = 7^ {asifx')) q{N, Z^') , (2) 

which leads to an exponential solution for the Mellin moments of the quark distri- 
bution, 



q{N, /i ) = exp 



.0 /^" 



■7jv (as(/i'^)) 



q{N,4)- (3) 



Note that here we express the solution in terms of a generic initial condition at 
some reference scale, as appropriate for the evolution of physical, measured parton 
distributions. When one instead considers parton-in-parton distributions, defined in 
QCD in terms of matrix elements of bifocal operators, one can use dimensional reg- 
ularization to express the solution as a pure exponential (with no prefactor), using 
the fact that the dimensionally regularized coupling vanishes with the scale [5], El] . 
Within the framework of dimensional regularization and in a minimal subtraction 
scheme, the structure of the anomalous dimension 7Af(as) at large values of A^ (cor- 
responding to the z -^ 1 limit) is known |22] to be single-logarithmic. It is of the 
form 



7;v («.) = -A (a,) InN + Bs (a.) - C^ (a.) ^r + ^^ («^) T? + ^{th] (4) 



InA^ . . 1 .^ / 1 \ 



where the function A{as) is one half of the cusp anomalous dimension 'yxi'^s), and 
N = Ne'^'^. The DMS proposal is that the functions C.y{as) and D^{as) are not 
genuinely independent, but they can be derived from the knowledge of A{as)- In 
turn, A{as) can be interpreted as a definition of the coupling in a suitable scheme, 
which has been variously described as 'physical', or 'bremsstrahlung', or 'Monte 
Carlo' scheme [23] ■ In order to implement this idea, DMS propose to replace Eq. ([T]) 
with 



/z^ A ^(^, /i^) = ^ ^ ^ (f , ^V) V [z, a. (^ 



(5) 



Here we have denoted by iplxjfi"^) a distribution which can be understood either 
as a fragmentation function or as a parton distribution; the parameter a = ±1 
serves to distinguish the two cases: cr = +1 for the space-like evolution of parton 
distributions, while a = —1 for the time- like evolution of fragmentation functions. 
DMS argue (and verify at two loops) that with Eq. ([5]) the evolution kernel is the 
same for both kinematics. Furthermore, at least up to second order in a^, the kernel 
V has no contributions at order (1 — z)°, so that it can be written as 

V {z, a,) = -jP^ + Bs (a,) 5(1 - z) + O {{1 - z)) . (6) 

(l-z)+ 



If one now chooses the cusp anomalous dimension (divided by the Casimir invariant 
of the appropriate representation, in this case Cp) as the definition of the couphng, 
setting y4(as(/i^)) = CFa'pn{^J^)1 one may conclude that all contributions to the 
evolution kernel that do not vanish as z ^ 1 appear at the first non-trivial order in 
this scheme. 

In the physical scheme, writing V{z, apn) = 'Pii.z) a^n/T^ + C^(«ph)' i^ i^ easy to 
construct an exponential solution, analogous to Eq. ([3]) but valid to NE order, for 
the distribution D. Indeed one may write 



/^ 






HN,fi' 



dzz''''Pi{z)apn{ — ]^ljiN,z''fi 



/^^ 



o-,.2n 



(7) 



The scale of the coupling can be shifted by using the (3 function, as 



2 d 



dij? 



^iN,fi' 



/' 

Jo 



dzz''-'Vi{z) 



«PH 

TT 



+ {l-z) Piapn) - a 



apH . 2 9 



TT 



/^ 



9/^2 



HN,fi'). 



One can now perform a Mellin transform, and introduce the anomalous dimensions 



^N~l 



7i(iV)= / dzz''~'P,iz) 
'o 



71 '(AT) = / dzz' 
'o 



~''-'[l-z)V^{z), (9) 

which clearly obey ^'{N) = ^{N) — ^{N + 1). One finds then 

.2\ 



^{N,f^' 



exp 



,,2 



Mo 



V^(iv,/i^) 



(10) 



l + a7'(iV) (apH(/x2)/7r) 
which is valid up to corrections vanishing as 2; — > 1. 

2.2 Moment integrals to 0{l/N) 

Let us now turn to the practical issue of evaluating the generic integrals appearing 
in the exponents of threshold resummations, to our required accuracy, i.e. including 
all correction of order 1/N . To this accuracy threshold-resummed partonic cross 
sections can be written as 



1 ^N-l _ 1 



Z)\. fill 



\n[a{N)\-H= dz— /,[ln(l-^)]+ / dz z^ -^ h[\n{l - 

Jo ^ ~ ^ Jo 

where H represents A^-independent terms. Expanding the functions f] in powers of 

their argument, as 

00 

/,[ln(l-z)]=^/i^)ln^(l-^), 

p=0 



(12) 



we can write 



In [aiN)] -H = J2 \fi'^ ^pW + f? J,{N) 



p=0 



(13) 



in terms of the basic integrals 



»1 N-l _ 1 

Vp(N) = I dz ^_^ lnP(l - z) , Jp{N) 



dzz^-^\nP{l~z). (14) 



In order to evaluate the integrals in Eq. ( IT4l) . we follow [20] and introduce two 
generating functions, defined by 



c.,..,.,,«--,a-,^-..mi-i 



(15) 



and by 



^^ ' ' Jo T{N + \ + l) N + \ 



[\Gv{\N) + l\, (16) 



lYioia. these definitions, one sees that the integrals in Eq. ( TT^ are given by 



V,{N) 



QV 



GviX,N) 



JpiN) 



A=0 



dP 



Gj{X,N) 



(17) 



A=0 



In order to evaluate the integrals explicitly to 1/A^ accuracy, we only need the first 
correction to Stirling's formula for the P-type integrals. 



T(z) = e-' z'-^/^ V2^ ( 1 + 



I 



12z 



1 + C 



1 



leading to 



Gv{X,N) 



A 



r(l + A) f^ A(l-A) ^ 



ATA 



2N 



while for the jT'-type integrals it suffices to take 

r(i + A) 



GjiX,N) 



Ni+x ■ 



(19) 



(20) 



We note in passing that, to 1/A^ accuracy, there is a simple relation between the JT" 
and the T> integrals; in fact 



MN)^-^VAN) + o(l, 



(21) 



which follows from an identical relation between the generating functions, 

GAX, N) = -^Gp(A, N) + (^) . (22) 

A useful way to evaluate both sets of integrals in the large N limit is to map them into 
simpler integrals, where the dependence on A^ has been moved from the integrand to 
the upper limit of integration. This technique is well known [21 EHj, and we extend 
it here to 1/A^ accuracy. Let the generating function of cutoff integrals be 

Gl{X,N)^I dzil-z)^-' = . (23) 

It is then easy to relate this function to the functions Gv and Gj. Expanding 
Eq. (TTOll in powers of A one finds 

Gp(A, N) = -G,(A, AT) + f^ ^ A^-i ^ , (24) 

fc=i 



where 

This can be rewritten as 



(25) 



A=0 



^-(^' ^) = E ^(-l)'"5(I^^-^^' ^) • (26) 



k=0 

Using Eq. (|22|) one then immediately finds 



«-(^.^')4s^'-i''5(SsT^'<'-'"'- <^'' 



Eqs. f l26l) and fl27l) can be used to evaluate directly the V and J' integrals to the 
desired accuracy, and indeed we will make use of this explicit evaluation in Section O 
One finds 



fc=0 

P 



■7. = ^ Er'"(i) (-inivr' + o(!^ 1 , (2S 



where T^''^ is the /c'th derivative of the Euler gamma function. On the other hand, 
one can use Eqs. (J26l) and (!27l) to directly relate the logarithm of the cross section 
to a cutoff integral of the same functions /i and /2 appearing in Eq. (ITTi) . This is 



useful when one needs to correctly account for running coupling effects to all orders, 
as done in [21 [20] . To the present accuracy one can write 

iV ^ A;! ^ ^ 5(lniV)^+iyo 1 - -^ 

We now move on to applying these tools to the concrete example of threshold re- 
summation for the Drell-Yan and DIS cross sections. 



3 An ansatz for next-to-eikonal logarithms 

In order to include NE effects in threshold resummation formulas we propose to 
modify the exponents in three ways. First of all, following DMS, we include sub- 
leading corrections in the argument of the running coupling. Second, we change 
the boundary of phase space accordingly. Third, and most relevant, we interpret 
the leading- logarithm function A{as) as arising from collinear evolution, and thus 
replace it with a NE generalization dictated by the DMS equation. This is done in 
the following way. While Eq. ([5]) cannot be diagonalized by means of a simple Mellin 
transform, it is however possible, as pointed out by DMS, to map the kernel V{z, as) 
in Eq. (^ back to the conventional evolution kernel, order by order in perturbation 
theory, if one explicitly performs the shifts in the arguments of Eq. ([5]) by the action 
of differential operators. Indeed, one may rewrite Eq. ([5]) as 

/^^^^(x,/z^) = £^e-"^^(^("^^^-'^^)^(^,/i^) P(.,a.(/.^)) , (30) 

where one should note that dependence on the coupling is only through the kernel 
V, while explicit scale dependence arises only in the distribution ip. Expanding the 
exponential and the kernel V in perturbation theory one is led to an equation which 
can be diagonalized order by order. When solved in this way, the DMS equation 
can be understood as a framework to generate classes of higher-order contributions 
to collinear anomalous dimensions using low-order information. In this spirit, we 
will write conventional resummation formulas, but we will generalize the collinear 
evolution function A{as) by including all terms that are generated by the DMS 
equation. As we will see, this will lead to slightly different implementations for 
space-like and time-like kinematics. Let us now consider our two examples in turn. 

3.1 The Drell-Yan cross section 

We first consider the Drell-Yan hard partonic cross section in the MS factorization 
scheme, denoted u}{N). We propose to generalize the exponentiation of threshold 

8 



corrections in the following way. 



In 



Q{N) 



^DY (as(Q')) + dz 



^N-l 



D 



a. 



2/02 



[i-^YQ 



+ 2 



{l-zfQ^/z 



Q2 






z,as{q'^ 



(31) 



where for simplicity we have set the factorization scale fi% = Q"^. Here and below 
we adopt the convention that the 'plus' prescription applies only to singular terms 
in the expansion of the relevant functions in powers of 1 — 2;. In other words, for 
a singular function f{z) with Laurent expansion f{z) = J2'^=^i fn (1 ~ -2)") and for 
any smooth function g{z), regular as z — > 1, we define 



dzg{z) f{z) = / 



dz 



9{z)-9{l) 



+ 



dzg{z) f{z) - 



/-I 



(32) 



In Eq. ( I3T!) . J?-dy(«s) is responsible for the exponentiaton of A^- independent terms, 
in accordance with [1]. It comprises purely virtual contributions given in terms the 
quark form factor, and real emission terms, which were denoted by Fj^g (a^) in |1]. 
The single- logarithm function D{as) can also be related to form factor data, and 
to the virtual part of the coUinear evolution kernel Bi{as), as was done in 
according to 



D{as) = ABs{as) - 2G{as) + (3{as) 



^/^' ^"^ 



(33) 



where G is constructed from single pole contributions to the quark form factor, 
as described in |4j. Finally, the DMS-improved space-like coUinear evolution kernel 
Ps{z, as) is given in perturbation theory by Ps{z, a^) = ^^1 Ps (z) (ag/vr)", where 



P^-\z) 



1-z 



^{n) ^ ^(„) j^^^ -z)+D 



7 



(34) 



Here A^^^ and C^ are the perturbative coefficients of the functions appearing in 



^(^) ; 



Eq. dl]), while D is related to the perturbative coefficients of D^{as) by the simple 



t{") 



shift D = D^ + y4*^"); this takes into account the explicit factor of z multiplying 
74(0;^) in Eq. flM|) . which in turn is responsible for the inclusion of NE terms in the 



D 



ordinary evolution kernel. In our normalization, A^^^ = Cp, Cy ^ 
at two loops 



.(1) 

■"1 



{1) 



0, while 



^(2) 



^(2) 



cl 



- - C(2) J CaCp - -UfCp 

--CACF+'UfCp. 



Cf 



^2 



(35) 



Notice in particular that the DMS procedure has brought into the resummation 
exponent abehan-hke terms proportional to Cp at two loops. As we will see, these 
terms do indeed find a match in the finite order expansion of u){N). The ansatz (|3T]) 
can be written in form of Eq. (ITTl) . and evaluated using the methods of Section [21 
In Section HI we will compare the perturbative expansion of Eq. (l3Ti) . with the 
coefficients given in Eq. fl35l) . to the exact results of [17]. In both cases, one may 
write the expansion 



UJ 



oo 
j=0 



2n 



2n-l 



J2 a„„ln'"Ar+ J2bn 



m=0 



m=0 






+ ^(l^)> (36) 



and then compare the expressions for the coefficients a„m and hnm arising from the 
resummation to the exact ones. 

3.2 DIS structure functions 

We consider next the resummation for the DIS structure function F2{N), in the MS 
factorization scheme. Phase space and kinematics in this case are somewhat more 
complicated, since one has to deal with the final state jet, which is approximately 
massless near threshold, as well as with initial state soft and coUinear radiation. We 
propose to generalize the conventional resummation formula as 



In 



A(iV) = ^Dis (as(Q')) + / dzz 





{i-z)Qyz 



Q 



,Ps 



z,as{q^ 



N-l 



1 



1-Z 



B 



a. 



+ 



{l-z)Q^/z 



dq' 



6P 



(l-2)2Q2/z ^ 



z)Q' 



z,as{q^ 



(37) 



Here, as above, J-'bis{c(s) is responsible for the exponentiation of A^- independent 
terms. The case of the DIS cross section in the MS factorization scheme was not 
explicitly treated in Ref. [1], but it is easy to work out the relevant contributions from 
the information collected there. Indeed, one can reconstruct the structure function 
F2{N) from the moment space ratio of the Drell-Yan cross section computed in the 

MS scheme to that computed in the DIS scheme, both given in |4j, as F2 (N) = 



^(Ms )(^jY^^^(Dis)(^jY^_ One then easily verifies that J-'-dis{c(s) comprises a virtual 
part, given by the finite terms in the modulus squared of the space-like quark form 
factor, plus a combination of real emission contributions, which can be written as 
(-^MS i'^s) — -^Dis(as)) /2 in the notation of [1]. The single-logarithm function B{as) 
can be associated with the evolution of the final state jet. It is interesting to note 
here that B{as) can also be expressed in terms of form factor data, plus virtual 
corrections to the collinear evolution kernel, plus a total derivative of lower order 
contributions, just like the function D{as) in Eq. ([33D- Indeed, one verifies that 



10 



existing results up to three loops are consistent with 

B(as) = BsM - G{as) + /5(«,) -^FB(a.) , (38) 

das 

with easily computed perturbative coefficients for the function F-Q^ag). Eq. fl38|) is 
in keeping with the general results of Ref. [25], where it was shown, at the ampli- 
tude level, that all IR and collinear singularities in massless gauge theories can be 
constructed from combinations of eikonal functions with the virtual collinear func- 
tion Bs{as), up to total derivatives with respect to the scale. Finally, we turn to 
the second line of Eq. (137|) . There, we have used the fact that the integration over 
the scale q^ has a range that can be split into two intervals, which correspond to 
different physical sources of radiation. Scales between the factorization scale Q^ 
and the soft scale (1 — z^Q"^ correspond to Drell-Yan-like initial state radiation, 
while scales between the soft scale and the jet scale, (1 — z)Q^, correspond to the 
evolution of the final state jet. Accordingly, in the ffist range we use the same 
space-like evolution kernel Ps(z, a^) that was employed in Eq. fl3Tl) . while in the sec- 
ond range we use the time-like fragmentation kernel Pt{z, as)- One may then define 
5P{z,as) = Pt{z,as) — Ps{z,as), and thus get to Eq. fl37j) . The function 6P{z,as) 
begins at two loops, where it is given by 



5P(2)(^) = --Cl {A\n{l - z) + 3) + 0{1 - z) . (39) 

Once again, using the methods of Section [2], we can expand both the resummed and 
the exact results for F2{N) in powers of logarithms of A^, and in inverse powers of 

A^, as 



oo 

%{N) = Y, ( 

1=0 



2» 2n—l In'™ AT 

m=0 m=0 



-Ol'^). ,40, 



We can then compare the resummed and exact values of the coefficients c„m and 
dnm, up to two and three loops, using the results of p^fTO]. 

4 Discussion 

We begin by checking the behavior of our ansatz at the one loop level. This is 
not trivial, since we have not added new coefficients in the exponent at one loop, 
and the only sources of 1/A^ terms are the expansions of the T>p integrals, and the 
simple modifications of phase space. Using the one loop results for the functions 
A{as) and D{as), we find that for the Drell-Yan cross section the one-loop exact 
result is recovered, including all corrections down to 0{1/N). Specifically, expanding 
Eq. ( I3T1) . we find bn = 2Cf and 6io = 0, which is exact. Note that 6io vanishes as a 
consequence of a cancellation between subleading terms in the expansion of the T>p 
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integrals and the modified phase space boundary. For DIS, including the one- loop 
value of the function B{as), we find that du = Cp/2 is correctly reproduced, while 
the non-logarithmic term at 0{1/N) is underestimated: Eq. (!37j) yields dio = Cp/S, 
while the exact result is dio = 21/8 Cp- We take this as evidence (to be reinforced 
below) that our treatment of phase space for the final state jet is sufficiently precise to 
reproduce single NE logarithms, but not enough to fix NE constants (of course at this 
level non-factorizing effects for the observable, leading to a failure of exponentiation, 
at least in the form of Eq. fl371) , may also be a source of the discrepancy) . 
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Table 1: Comparison of exact and resummed 2-loop coefficients for the Drell-Yan 
cross section. For each color structure, the left column contains the exact results, 
the right column contains the prediction from resummation. 

At the two-loop level, we proceed as follows. Since our aim is to verify our ability 
to reproduce NE terms, suppressed by a power of A^, we include in the exponent all 
terms that are required to reproduce ordinary Sudakov logarithms, i.e. the two-loop 
values of the functions A{as) and D{as) for the Drell-Yan cross section, and of the 
function B{as) for DIS. We include the two- loop DMS-induced contributions C\ 

D and 6P^'^\z) as well, since they are responsible for effects that originate at two 
loops, and can only be reproduced by their inclusion. Our results are summarized 
in Tables 1 (for the Drell-Yan cross section) and in Table 2 (for the DIS structure 
function) . 
We observe the following. 

• The leading non- vanishing NE logarithms (In N/N for the 'abelian' terms 
proportional to Cp, and In N /N for non-abelian terms) are correctly repro- 
duced by the exponentiation, both for DY and for DIS, and separately for each 
color structure. 

• Next-to- leading NE logarithms (In N/N for terms proportional to Cp, and 
In N/N for non-abelian terms) are reproduced with remarkable accuracy for 
the Drell-Yan process (in fact exactly for the njCp color structure), and rea- 
sonably well for the DIS process. 

• The remaining NE logarithms, i.e. single logarithmic terms proportional to 
Cp, are well reproduced by exponentiation for the Drell-Yan process, but only 
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roughly approximated for DIS. Non-logarithmic NE corrections are not well 
approximated by the exponentiation. 

More specifically, we note that for the Drell-Yan process the only source of 
terms proportional to Cp In^ N/N is the DMS-induced coefficient C^ ; indeed, 
the fact that 610 = ensures that no such term can arise from the square of 
the one-loop contribution. This contribution, yielding 622 = 4, is an excellent 
approximation to the exact result, 622 = 7/2. For DIS, as might be expected, 
the situation is somewhat more intricate; indeed ^22 receives contributions 
from three sources: the square of the one-loop exponent, C\ , and 5P^'^\z); 
also here, however, the final result, ^22 
exact answer, (^22 = 39/16. 



55/16, is a fair approximation of the 
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Table 2: Comparison of exact and resummed 2-loop coefficients for the DIS structure 
function. For each color structure, the left column contains the exact results, the 
right column contains the prediction from resummation. 

Clearly, since some of the DMS modifications enter the stage at two-loops, our re- 
sults verify that these contributions improve the approximation, but do not really 
test exponentiation. We can put at least our DIS ansatz to a more stringent test by 
comparing to the complete three-loop calculation performed by Moch, Vermaseren 
and Vogt [12]. In this case, since our aim is to test exponentiation at NE level, 
we have included the three-loop value of the function B^ag)-, contributing to single 
Sudakov logarithms, but we have not included three-loop DMS-induced contribu- 
tions such as C\ and 5P^'^\z). We can then expect reasonable agreement only for 
a limited set of NE logarithms. Since at three loops one finds six independent color 
structures, up to five powers of NE logarithms, and transcendentals up to C5, we 
do not include here the lengthy tables of coefficients, but we give the most relevant 
results. 

The three-loop analysis confirms that leading non-vanishing NE logarithms (in 



this case In N /N for the color structure C|., In N /N for the color structures CaC^ 
and rifCp, and In^ N/N for the color structures Cj^Cp, ri'lCp and UfCACp) are 
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exactly reproduced by our resummation ansatz. Next-to-leading NE logarithms are 
reasonably well reproduced: specifically, for all color structures and separately for 
each degree of transcendentality the approximate results from the resummation have 
the same sign and similar numerical values to the corresponding exact results. In 
particular, this applies to the coefficient d^, whose exact value is 57/64, while the 
approximate result is 109/64. Since ^34 arises in part from interference between the 
NE coefficient C) and the leading one-loop Sudakov logarithms in the exponent, we 
take this as mild evidence in favor of the exponentiation of DMS-induced corrections. 
To summarize, we have provided an ansatz to include in threshold resummation 
a set of next-to-eikonal corrections, allowing for subleading phase-space effects, and 
including the modified collinear evolution proposed by Dokshitzer, Marchesini and 
Salam. It is understood that these modifications of conventional threshold resum- 
mation do not exhaust all possible sources of NE threshold logarithms, and indeed it 
may be expected that some such corrections might break Sudakov factorization and 
fail to exponentiate. By comparing our ansatz to finite order perturbative results for 
the Drell-Yan and DIS cross sections, up to three loops, we have however provided 
evidence that at least the leading non- vanishing NE logarithms do indeed exponenti- 
ate according to our proposal. We have furthermore provided evidence that the DMS 
equation induces a definite improvement for resummation at NE level: for example, 
abelian-like next-to-leading NE terms that conventional resummation completely 
fails to generate are accurately approximated when DMS evolution is implemented. 
In general, it is clear that our ansatz gives better results for the Drell-Yan pro- 
cess, presumably thanks to its simple phase space and kinematics. The presence of 
the final state jet in DIS, and the related constraints on phase space, may require 
a more detailed factorization analysis in order to collect all sources of NE terms, 
and indeed may well induce a breakdown of simple Sudakov factorization at NE 
level. To aid this preliminary exploration of NE exponentiation, we have provided 
here some practical tools that will be useful in future extensions of this work, and 
we have taken the opportunity to note a connection, given in Eq. (l38l) . between 
the jet function B{as) and the virtual collinear function Bs{as), as was previously 
done for the soft function D{as) in the Drell-Yan cross section [23]. We believe 
that this work provides further motivation both to include leading NE correction 
in phenomenological resummation studies, and to pursue the corresponding theo- 
retical work. Indeed, a full understanding of NE threshold logarithms must await a 
thorough analysis of soft gluon radiation beyond the eikonal approximation in the 
non-abelian theory, and specifically an adequate implementation of Low's theorem, 
mapping its boundaries of applicability in the case of massless QCD. 
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